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In this study, the flow characteristics of air around a 
circular cylinder were determined by two means: by 
analytical solution of the Na vier-Stokes equations, and 
actual measurement of the flow itself. A DISA LD 5586 laser 
doppler anemometer system was employed to measure normal {v) 
and parallel (u) velocity components of a laminar flow 
regime about a three-sixteenths inch diameter cylinder 
mounted in a 32 x 45 inch wind tunnel (Reynolds number 
approximately 40) . 

A finite element grid was constructed using triangular 
elements which encompassed the upper half-cylinder. Flow 
measurements using the laser velocimeter were confined to 
this region cf interest. 

Theory was compared against experimental results, and 
the feasibility of the utilization of this particular type 
of laser anemometer was evaluated. 
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II. 



DESCRIPTION OF LASER ANEMOMETER 



A basic laser anemometer is comprised of four primary 
elements: 



(1) a coherent light source 

(2) optics 

(3) signal converter (photomultiplier tube) 

(4) signal processor/display 

The principle of the laser anemometer is to create a 
grid pattern and measure the speed with which a particle in 
a fluid free stream passes through the grid. The following 
assumptions were taken into consideration: 

(1) particles injected into the free stream are 

small enough to faithfully follow stream 
lines of the flow, thus by measuring particle 
characteristics, we are actually measuring 
the medium's characteristics. This is easily 
achieved by use of particles which are on the 
order of a few microns in diameter. 

(2) particle generation results in a consistent 

particle size, or at least a majority of the 
particles are within a specified range of 
sizes. For example, DOP (di-2-ethylhexyl 
phthalate) particle generators produce 

particles with a mean diameter of 0.75 
microns and burning turbine oil produces 
particles on the order of 4-6 microns (from 
electron microscope) . 
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A laser is incorporated because of its coherency and 
power. The aforementioned grid is in reality the result of 
two coherent beams interfering with themselves (Fig. 1) . 
The optics produce the beam pair by conventional splitting 
and focusing. Because the laser beam intensity is of a 
Gaussian distribution the grid created by constructive and 
destructive interference is not uniform over the entire 
focal volume. The fringes created on the outer portions of 
the horizontal axis are not perfectly formed because of 
slightly different intensities in that area. The main 
portion of the interference grid is uniform and, most 
importantly, the fringe spacing (D^) is determined only by 

the laser wavelength (lambda) and the half angle of the beam 

pair convergence (theta/2) . 



= lambda/ (2*sin (Theta/2) ) 



The number of fringes (N^) may also be computed, 
provided the beam focusing convergence angle (delta theta) 
is known ,by the relation; 



N^=8*tan (Theta/2) / (PI*Delta Theta) 



The number of fringes is important only to ensure that the 
signal processing apparatus will receive enough information 
to . make possible further data reduction. If an insufficient 
number of fringes exists, not enough signal will be 
available for velocity determination. 

Solid particles with a mean diameter which is less than 
half the fringe spacing are introduced into the free stream 
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Figure 1 
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of the fluid to he measured. When the particle is small 
compared to the fringe spacing the scattered light exhibits 
essentially 100% modulation, while particles equal to the 
fringe spacing produce negligible modulation[ Ref . 2]. The 
particles may be produced by atomization of water, dop-air 
reaction, burning oil, smoke bombs, or may be bought 
commercially, as e.g, polystyrene spheres and some solid 
chemicals such as zinc oxide. Choice of the particle to be 
used is normally dictated by the fringe spacing in the focal 
volume, and the nature of the medium to be measured. 
Regardless of choice, the particle must be small enough to 
follow the streamlines as closely as possible, because the 
laser anemometer actually measures the particle velocity and 
unless the particle's motion is a true representation of the 
fluid velocity, pertinent data may not be obtained. 

Inside the measuring (or focal) volume, bands of 
darkness and brightness exist due to the interference of the 
two beams. Particles entering the focal volume will 
alternately scatter light as they pass through regions of 
constructive interference ( typical focal volume dimensions 
are 300 microns wide by 4 millemeters long and contain 
approximately 64 fringes) . A photomultiplier tube located 
on the optical axis and focused into the measuring volume 
receives this scattered light and converts what appears to 
be a modulated light signal into a series of voltage pulses. 
Each voltage pulse correlates to a particle crossing one 
fringe. Eecau-se the spacing of the fringes is known, the 
frequency of the voltage pulse train may be equated to the 
velocity of the particle by the following relationship: 



f=V /D = 2*V *sin (Theta/2) /lambda 
of o 

where D is the fringe spacing and V is the velocity of the 
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particle. Now the processor is left with a fairly simple 
task of counting the frequency and multiplying it by a 
calibration factor, which is, in essence, a fringe spacing 
parameter to arrive at a velocity. 

As long as we are assured that only one particle is in 
the measuring volume at one time, the apparatus will operate 
correctly, but this is an unrealistic assumption. In order 
to achieve such a situation, the seeding density of the 
particles would be exceptionally sparse and the resulting 
data rate would be so slow as to preclude timely data 
reduction. To overcome this obstacle, a validation circuit 
is incorporated into the counter/processor. When a particle 
crosses a fringe spacing, two counters are triggered. One 
measures the time duration for five fringe crossings and 
holds this information. The second counter continues to 
measure the time for a total of eight fringe crossings. The 
time of the first counter is then compared with that of the 
second counter, and if the first time is indeed five-eights 
of the second time (within specified limits, i.e., comp, 
accuracy) a validation gate is opened and the resultant 
frequency is then stored into an ensemble cell. The ensemble 
cell is adjustable to hold various numbers of counts an^ 
average them for the final velocity display. The final 
displayed velocity then is that component of fluid velocity 
which is perpendicular to the optical axis and contained in 
the same plane as the beam pair (Fig. 2). 

The anemometer may be operated in three modes: 
reference beam, forward scatter, and back scatter. The last 
two modes are most frequently employed and are described 
here: 

Forward scatter mode: in this mode, the photomultiplier 
tube is located on the optical axis and focused in a 
direction directly opposed to the direction of the emerging 
beam pairs. This mode offers the strongest attainable 
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signal. 

Back scatter mode: again the photomultiplier tube is 
mounted on the optical axis, but now it is aligned in the 
direction of the emerging beam pairs. This mode affords the 
ability tc measure fluid velocities in confined spaces or in 
situations where the dectector cannot be mounted in forward 
scatter, such as in measuring velocities in a compressor, 
additionally, operation in this mode requires a higher 
powered laser due to the reduced intensity of back 
scattered light (Fig. 3) . 
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ANALYTICAL DERIVATION 



The finite element method is a numerical analysis 
technique fcr obtaining approximate solutions to a wide 
variety of enginering problems. It offers a means to solve 
complex continuum problems by sub-dividing them into simpler 
interrelated problems. The nature of this engineering 
problem does not readily lend itself to a closed form 
solution. Therefore a numerical solution is dictated. 

The method of finite elements is well employed where 
there are irregular geometries or unusual specifications of 
boundary conditions. The method gives a piecewice 
approximation to the governing equations. Figure 4 shows 
the finite element grid employed in this flow problem. 
Approximating functions (interpolation functions) are 
defined in terms of the values of the field variables at the 
nodes.'^For the finite element representation, the nodal 
values of the field variables become the new unknowns and, 
once found, the interpolation functions define the field 
variable throughout the field of elements. The degree of the 
interpolation polynomials is governed by the number of nodes 
per element and dictated by the order of the governing 
equations. 

The basic principle of the finite element method is to 
divide the solution domain into a finite number of 
subdomains (elements) which are connected only at node 
points. Interpolation polynomials are created for each 
element and a collecton of interpolation functions over the 
whole soluticn domain provide a piecewise approximation to 
the field variable. 

Introduction of generalizations facilitate direct 
derivation of finite element equations from governing 
differential equations. The method of weighted residuals 
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(Gaierkin's Kethod) is a technique for obtaining approximate 
solutions to non-linear partial differential equations. We 
first assume the general function behavior of the dependent 
field variable in some way so as to approximately satisfy 
the given differential equations and boundary conditions. 
After applying this approximation to the original equations, 
seme error will result, which we desire to minimize or drive 
to zero over the entire solution domain. 

The basic set of differential equations governing the 

problem from a velocity and pressure formulation approach 

are (after linearization by approximating non-linear viscous 

terms u , v ) : 
n n 



dv _ 

continuity ^ = O 



^ II ^ V/ J. = 

momentum Qy “ ”|0 



a 



0x ay 



u , V , and E are approximate solutions to the flow problem 
n n n 

th 

and may be visualized as the n iteration of the procedure 

of solution. These equations used in formulation of the 
problem are by no means unique. Stream functions and stream 
function with vorticity formulations have also been employed 
to establish finite element solution procedures. 

Inasmuch as the velocity terms are of second order, a 
second order interpolation polynomial must be used to 
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calculate these unknowns. The pressure terms are of first 
order and hence require only first order interpolation. 



The basic element chosen for this solution was the 
triangle and, in a localized scheme, is visualized as 
follows : 




Nodes 2, 4 and 6 lie equidistant between their respective 

corner nodes. 



Nodes 1, 3, and 5 are referred to as corner nodes and 

are the basis for the linear interpolation polynomials 
derived fcr the pressure terms. For any triangular element, 
the corner nodes are described by their Cartesian 

coordinates x , y . 

i i 

Consider the following local element: 







The coordinates of node one are x , v ; node two x , y ; 

11 2 2 

etc. and the interior point (x,y) may be described as a 
linear combination of the three corner nodes: 
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x=L X +L X +L X 

1 i 2 2 3 3 



y=L y +L y +L y 

11 2 2 3 



3 



wb^re L is the natural coordinate (interpolation 
i 

polynomial, here linear) or weighting function which relates 
the coordinates of the corner nodes to any interior node. 
L has a value of one at the point x ,y and is zero at the 

i ^11 

other two nodes of the element. 



The L 's expressed as linear functions of the Cartesian 
i 



coordinates then are 



\ 



Lj(x,y)= (a^+b^x+c^y) 

tSt (a +b x+c y) 
2 2 2 2 



where 2.^ = 



1 X y 

11 

1 X y 

2 2 

1 X y 

3 3 



which is equal to the area of the triangular element and 
a =x y -X Y , b =y -y , and c =y -y , the determinants of 

1 2 3 3 2 12 3 112 

the minors corresponding coefficients a^, b^, or c^, and the 



other coefficients are similarly computed. 



For the velocity terms, second order interpolation 
polynomials are indicated, while the pressure terms remain 
first order. Hence the introduction of more nodes (6) will 
result in a new set of weighting functions which may also be 
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interpolation 



expressed in terms of the natural coordinate 
polynomials; 



N =2L -L 

1 11 

N =2L^-L 

3 2 2 

N =2L^-L 

5 3 3 

P 

N =L 

i i 

e 4 

u = X N u 
• i i 



N =4L L 
2 12 

N =4L L 

4 2 3 

N =4L L 

6 13 



e 4 

V = 2 N V 

» i i 



e a P 

P = P (xxxx) 
» i 



Application of Galerkin's Method to the governing 
equations results in the following intergral expressions; 



L 



N-r(lrx‘iy) dB. = o 



here JX 



implies "over the area of the triangle". 



Ihen we may make the substitution of the field variable 

in terms of x and y into terms of L . 

i 



e.g. ; 


gxven 


then 


au _ 




dx “ 


where 


3Li 




ax ■ 


and a 


similar 



u(x,y) <=> 
3L; 3x 

hi 

2A 

procedure for 



u(L.) 

X 




te) 3^ 9 Li 
^ dL\ 3y 
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where 



3U . C; 

" Z.A 

e e e 

Substitution of equations (xxxx) into u , v and P and 

application of the above chain rule allows integration by 
parts : 



K = 

1 



>4 / . 

ax 






3y Sy 





y 





The evaluated integral expressions then form the influence 
matrix as fellows: 



K ! 0 1 


K 


r 

1 

1 

T 

1 


a 


0 ! K * 


K 


l_ _ll 


_ 3 


K ' K 1 


0 


l_ 2 1 3, 





Specification of the field variables or right hand sides of 
the matrix equation then provides a defined problem: 



Additionally, evaluation of the integral expressions yields 
boundary terms which must be computed for nodes which lie on 
the boundary of the control surface: 



i = \ N 
1 ) i 



* 

I ds 



R 






♦ 

I ds 






* * A * <jeJ ^ Q 

where X and I =^^W*7l-ny.r 

It should be noted that along vertical boundaries there is 
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no outward pointing normal n and similarly along horizontal 

y 

boundaries no n exists, and furthermore these terms are not 

X 

computed where u and v are specified. Additionally, it was 
assumed that along the left-and right-hand side and upper 
boundaries and. ^/^'j were identically equal to zero 

because these boundaries were placed at five • or more 
characteristic lengths (diameters) from the cylinder and 
considered to be outside the realm of influence of the flow 
about the cylinder. 

The terms comprising K^, and are entered element 

by element into the main program in steps FHJ02080 through 

QZ00170. Then each element's matrix is compiled into the 
field matrix (comprised of all elemental matricies) and is 
solved iteratively to output the horizontal velocity 
component u, the vertical component v, and the pressure P at 
every node in the field. 
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IV. 



D ESCR I PTION OF EXPER IMENT 

& right-angled, three-sixteenths inch cylinder made of 
stainless steel was mounted in an AEROLAB series 90 32x45 
inch low speed wind tunnel through the floor of the test 
section (Fig. 5) . The portion of the rod protruding down 
through the test section was mounted on a 
two-degree-of-freedom, micrometer controlled bed which 
allowed the rod to be maneuvered in the vertical plane 
parallel to the flow direction. This scheme enabled the 
laser velocimeter and photomultiplier tube to be mounted, 
aligned and left undisturbed for the entire period during 
which data were taken. The bed was positioned to read 
locations in the control volume which directly coincided 
with the location of the nodes in the finite element mesh 
(Fig. 6) . 

In the interest of precluding as many opportunities for 
error as possible, all u velocity components were first read 
and then components in a forty-five degree position were 
taken thus necessitating only one re-orientation of the beam 
pairs. 

The plexiglass window prefacing the optics transducer 
of the test section had one-quarter inch holes drilled into 
it to allow the beam pairs to pass unimpeded, thereby 
minimizing distortion of the focal volume fringe pattern. 
The photomultiplier tube was likewise accommodated through 
the far window. 

The seeding mechanism employed zinc-hexaclorcthane 
based (non-toxic) white smoke bombs. The nature of the 
closed return wind tunnel permitted continuous measurement 
without interruption for re-seeding of the flow. By means of 
electron microscopy, the average size of the smoke particles 
was determined to be on the order of two microns (Fig. 7) . 
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Cylinder mounted in wind tunnel and LDV setup. 



Figure 5 
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Traversing mechanism 



Figure 6 
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Zinc-hexachlorathane particles 
(X5000 lcm= 2.^m) 



Figrire 7 
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The final velocity component measured at each node was 
the result of averaging three sucessive readings from the 
velocimeter processor. Each of these readings represented 
the average of 4096 readings which had been validated to 
within 1.5% of the exact five-eigh.ths validation count 
criteria. 

The flow in the wind tunnel was stabilized at five 
inches per second to maintain a Reynolds number of about 40, 
which is in the realm of creeping flow. The experimental 
setup was designed to measure the steady flow about the 
cylinder, hence the relatively low Reynolds number, and no 
effort was made to obtain unsteady flow measurements. Three 
separate runs were conducted, two over the forward half of 
the control volume and one over the rear half. Because the 
beams of the laser anemometer were approximately 
one-sixteenth inch in diameter, several node position 
measurements were excluded due to inaccessibility. The 
majority of these were in close proximity to the cylinder. 
In order to obtain both an x- and a y-velocity component of 
the flow at every node which was accessible, two 
measurements were made, but because of the very small 
vertical component present in most of the flow, direct 
vertical beam pair measurement was unsuitable. Instead, 
after obtaining the horizontal component for each node in a 
run, a second series of measurements was made with the beam 
pairs inclined forty-five degrees from the vertical. In this 
fashion the y-component was determined by the following 
relationship: 
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V 










• N 
I N 



It can fce seen here that the unknown resultant (K) has 
a projection onto each axis representing the orientation of 
the beam pairs, thus the resultant is the vector which 
emanates from the origin and terminates at the intersection 
of the perpendiculars constructed from the u and v(at 45 
degrees) vectors. 
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V, 



D BSCHIPTION OF COMPDTER PROGRA M 



The program found in the appendix has several ccmment 
cards inserted for clarification. However, a broader 
presentation is offered here. The program may be subdivided 
into four main sections: 

Section 1: dimensioning and zeroing of 
parameters and arrays. (FLD00020 through 
FLU00690) 

Section 2; input of node coordinates, system 
topology and known values of u, v, P, and 
right hand sides, and output of the same. 
(FL000700 through FLU02090) 

Section 3: computation of interpolation 
polynomials on an elemental level and 
arrangement into local element arrays. 
(FLU021C0 through GO00320) 

Section 4; loading of the global influence 
matrix (TM) and solution of (TM) (T) = (Q) by 
Gauss-Jordan partial pivoting and back 
substitution, with an iterative routine 
incorporated until convergence is attained, 
and associated print out. (FLU04470 through 
FLU04880) 

The computer routine systematically tested the solution 
of field variables with the previous solution until all 
terms in the T and T1 arrays agreed within 0.001 (ESPILA) . 
If the convergence test failed, the process was repeated 
from section 3. 

Initially, the program was run on an IBM 360 in FORTRAN 
CLG mode, then an ensuing run was completed on the FORTRAN 
HCLG compiler and a time reduction of two-thirds was 
realized. The program, as given in the appendix, required 
approximately 375K bytes of core to compile and 675K bytes 
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t. 












to execute and output. Substantial space savings may be 
realized by utilizing banded matrix storage methods vice the 
method employed here, but due to time limitations this 
avenue was not fully explored. 

The program attempts to solve simultaneously for all 
the field variables. However, introduction of zeros on the 
main diagonal of the TM influence matrix resulted in a 
singular matrix. Another approach was explored where the 
equation of continuity 

|y = O 

was replaced by 

, f §y 3u _3u9yl 
2 1^ Qx 3y 3y J 

and the resulting terms were added to the TM matrix and 
their respective right hand sides. The only field variables 
affected by this form of the continuity equation were the 
pressure variables, and depending on whether the lacal 
corner node involved in the local pressure equation was on 
the boundary or interior, different right hand sides had to 
be evaluated. This is the purpose of steps EL004460 through 
G000320. The simultaneous solution technique for all field 



variables applied 


to this 


new set of 


equations did not 


converge. 










A third, unimplemented technique 


may be considered. 


Given the original 


matrix 


equations: 




(1) 


K 0 


+ 0 


+ K P = 


R 


1 




2 


1 


(2) 


0 


+K V 


II 


R 




1 


3 


2 


(3) 


K 0 


+ K V 


+ 0 = 


0 


2 


3 






Solution for u and 


V will 


result 


in 








-1 


- 1 




(1) 


0=-K 


K P 
1 2 


+ K R 

1 1 
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m 




II 

w 






■% 




( 2 ) 



V=-K K P +K 



1 3 



1 



R 



2 



and multiplication by and yields 



3 



K 1>=-K K K P + 
2 2 1 2 



K K R 



2 1 1 



K 7=-K K K P + 



K K R 



3 3 1 3 



3 1 2 



Summation of these equations allows substitution for the 
third matrix equation; 



The matrix operating on P may be calculated on an 

algebraic level and entered, in terms of a , b and c , into 

i i i 

the program in order to solve for the pressures at every 
corner node. The calculated P's then would be entered into 
the field variable array and the partial pivoting routine, 
which involves only those terms where the field variable is 
not specified, would operate primarily on the first two 
matrix equations. In previous computer trials where all 
corner node pressures were specified, although not known to 
be correct, calculations were properly performed. This 
method would allow the actual pressures present in the 
control surface to be entered into the computer routine for 
further calculation of the velocity components. 




+ K K R 



3 1 2 
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VI. ANALYSIS OF RESTJLTS 



No computer routine output was obtained for comparison 
of results, due to the inability of the program to converge. 
As such, comparison of the two methods cannot be offered 
here. The non-convergence of the program was investigated 
from several aspects: validity of the non-linear terms, 
correctness of the loading of the TM matrix on both the 
elemental and global levels, correctness of input, 
performance of the partial pivoting subroutines and the 
iteration scheme; all without success. In test cases the 
program initially solved for the unknowns within magnitude 
accuracy for a majority of terms. However, sucessive 
iterations resulted in oscillating solutions, even after 
twenty iterations. Convergence was expected after ten to 
fifteen iterations without under- or overrelaxation factors 
incorporated into the iteration scheme. After the new 
formulation of the continuity equation was incorporated into 
the computer routine, the program execution again resulted 
in a singular matrix when only the left and right hand and 
upper control surface boundaries had their respective 
pressures specified. Underflow occurred in the execution of 
the subroutine SOLVE at step SOL00260. Subsequent 

th th 

investigation revealed that the 178 and 179 rows of the 
reduced matrix equations were the cause of the underflow. 

The data which were taken for runs number two and three 
are displayed in the appendix and are plotted on the finite 
element grid in Figure 8. Without the analytical results to 
compare the measured data only a qualitative statement may 
be made pertaining to the LDV output. It is readily seen 
that the trend of flow directions around the cylinder 
appears to be correct. By and large, the flow in close 
proximity to the cylinder was accelerated over the top and 
decelerated near the stagnation points. It is pointed out 
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ii 









here that the plot is a composite of two different trials, 
and as such, different free stream velocities were present 
for each run. In Figure 8 the vector representing the flow 
at a node was plotted and the magnitude of the vector in 
inches per second was displayed above the vector or by the 
node from which the velocity vector emanated. 
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Velocity data plut 
Figure 8 




k. EHEOB ANALYSIS 

ae-positioning of the cylinder to its base location was 
accomplished with relatively low position error. Prior to 
each run the cylinder was aligned using the beam pairs as 
reference. In the horizontal mode, the beam pairs would 
strike the leading and trailing portions of the cylinder and 
reflect to positions horizontally opposed on the far side of 
the wind tunnel ^when correctly zeroed vertically. The same 
procedure was repeated for the vertical beam pair except 
position was checked by the upper beam reflecting to a 
vertical position on the far wall and the lower beam (at 
reduced intensity) striking the cylinder support on 
centerline. This procedure was repeated several times and 
the maximum deviation between any aligned positions was 0.2 
millimeters horizontally or vertically. Thus, ir was 
determined that the maximum positioning error inherent in 
the system would be the ratio of the smallest coordinate of 
any node measured and the aforementioned 0.2 millimeters or 

0.2/3.75=0,5% 

Preestream velocity measurements were likewise tested 
for fluctuations. Typical measurement variations ranged 
from 0.129 to 0.142 (run number three; freestream 
measurements) in the horizontal mode and 0.0947 to 0.100 in 
the forty- five degree mode (velocities in meters/seconds) . 
These figures represent a total of 10% variation which may 
be considered as ± 5.% from the mean. 
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VII . CO NCLUSION 

The experiment, as described, reguireTi approximately 
two to four hours for each run. aside from the initial 
alignment of the components and the cylinder, this was 
perhaps the greatest drawback of the system. Next, the 
exclusion of regions in close proximity to the cylinder from 
measurement was deemed to bfe a major factor in the small 
amount of data actually taken. The latter may be improved 
with a different experimental arrangement; particularly, the 
plexiglass section directly prefacing the optics transducer 
should be completely removed, thus allowing all angles of 
beam pair inclination, and in turn availing more locations 
to measurement. Another improvement which could be 
implemented would be unitization of the laser, transducer 
and photomultiplier tube arrangement. This would decrease 
the LDV's sensitivity to vibration and provide a more 
accurate optical axis alignment. The aforementioned 
improvements would still not decrease the time reguired for 
data collection, but the introduction of a two component 
system [Ref. 4] would allow simultaneous measurement using 
orthogonal beam pairs. Furthermore, this type of system 
would not concern itself with the errors introduced in beam 
pair re-alignments performed with the single beam pair LDV. 

The reproducibility of the experiment was judged to be 
good with the exception of exact duplication of the wind 
tunnel Velocity between successive trials. This, of course, 
is primarily a function of the controllability of the tunnel 
motor and propellor assembly. It is estimated that with 
adeguate facilities for incorporation of some of the 
aforementioned improvements, more accurate data may be 
obtained. The LDV displays great potential as a means of 
measuring a wide variety of dynamic fluid regimes accurately 
and without introduction of disturbance to the flow itself. 
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APPENDIX A 



LDV DATA 



RUN #2 



Node 


U(cm/ sec) 


\^45°(cm/sec) 


R( in/sec) 






12.43 


8.43 


4.90 


^3 


65 


12.33 


8.49 


4.85 


-1.5 


66 


13.10 


8.33 


5.18 


-5.8 


67 


13.17 


8.25 


5.22 


-6.5 


68 


13.87 


8.72 


5.49 


-6.3 


69 


13.87 


8.56 


5.50 


-7.2 


73 


14.03 


8.58 


5.57 


-7.7 


75 


13.40 


8.83 


5.29 


-3.9 


76 


13.43 


8.62 


5.31 


-5.3 


77 


13.00 


8.62 


5.13 


-3.6 


78 


13.47 


8.87 


5.32 


-4.5 


79 


12.87 


8.58 


5.08 


-3.3 


80 


13.63 


8.91 


5.38 


-4.3 


81 


13.20 


8.02 


5.25 


-8.0 


96 


10.50 


6.94 


4.14 


-3.7 


97 


12.83 


8.64 


5.06 


-2.7 


98 


13.47 


9.54 


5.30 


0.1 


99 


13.27 


8.92 


5.23 


-2.8 


100 


13.43 


9.25 


5.29 


-1.5 


101 


12.33 


8.83 


4.85 


0.7 


102 


9.77 


6.69 


3.85 


-1.8 


105 


13.03 


9.32 


5.12 


0.8 


106 


11.93 


8.68 


4.70 


1.7 
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1 

2 

3 

4 

6 

7 

8 

9 

13 

14 

15 

16 

19 

20 

21 

22 

23 

24 

25 

26 

29 

30 

31 

35 

36 

37 

38 

39 

54 

55 

56 

57 

58 

59 



U( cm/sec) ^@45°(cm/sec) R( in /sec) 

5.44 






^ 9C. 



13.73 


10.77 


13.70 


10.43 


13.77 


10.33 


13.23 


10.03 


13.77 


10.93 


13.80 


10.40 


13.87 


10.17 


13.27 


10.08 


13.50 


10.47 


13.57 


10.40 


13.53 


10.20 


12.30 


9.93 


13.70 


10.43 


14.60 


10.87 


13.83 


10.57 


13.60 


10.06 


13.57 


10.23 


12.93 


10.40 


12.50 


9.92 


11.90 


9.64 


13.13 


10.16 


12.57 


9.80 


11.30 


9.76 


14.17 


10.47 


14.03 


10.43 


13.23 


10,07 


12.97 


10.09 


11.83 


10.03 


14.03 


10.87 


14.17 


11.03 


15.20 


10.57 


15.17 


10.47 


16.27 


10.00 


14.37 


10.20 



5.43 


3.49 


5.22 


4.11 


5.46 


7.00 


5.44 


3.77 


5.46 


2.11 


5.24 


4.27 


5.34 


5.5^ 


5.36 


4.80 


5.34 


3.76 


4.89 


8.05 


5.41 


4.38 


5.76 


3.02 


5.46 


4.63 


5.36 


2.65 


5.35 


3.79 


5.14 


7.84 


4.96 


■ 7.00 


4.74 


8.27 


5.19 


5.40 


4.98 


5.86 


4.56 


12.50 


5.58 


2.60 


5.53 


4.36 


5.22 


4.36 


5.13 


5.72 


4.75 


11.24 


5.55 


5.58 


5.61 


5.76 


5.98 


-0.90 


5.97 


-1.36 


6,46 


-7.46 


5.66 


0.20 
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